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Abstract. For a cylindrical superconductor surrounded by a normal mate- 
rial, we discuss transition to the normal phase of stable, locally stable and crit- 
ical configurations. Associated with those phase transitions, we define critical 
magnetic fields and we provide a sufficient condition for which those critical 
fields coincide. In particular, when the conductivity ratio of the supercon- 
ducting and the normal material is large, we show that the aforementioned 
critical magnetic fields coincide, thereby proving that the transition to the 
normal phase is sharp. One key-ingredient in the paper is the analysis of 
an elliptic boundary value problem involving 'transmission' boundary condi- 
tions. Another key-ingredient involves a monotonicity result (with respect to 
the magnetic field strength) of the first eigenvalue of a magnetic Schrodinger 
operator with discontinuous coefficients. 
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1. Introduction 

A type II superconductor placed in an applied magnetic field undergoes three 
phase transitions: When the intensity of the applied field is below a first critical 
value ifci , the sample exhibits the Meissner effect and remains in a superconducting 
phase. When the field is increased further past Hc^ , the sample is in a mixed state 
and the magnetic field penetrates the material in quantized vortices. Increasing 



2000 Mathematics Subject Classification. Primary 35Q55, Secondary 35J25, 35A15, 58E50, 
81Q10. 

Key words and phrases. Ginzburg-Landau type equations, Schrodinger operator with magnetic 
field, semiclassical analysis. 

1 



2 



S. FOURNAIS AND A. KACHMAR 



the field further past another critical value , the sample exhibits surface super- 
conductivity, and when the field is increased above Hcs, superconductivity is lost 
and the sample is in the normal phase. The above picture is rigorously established 
for extreme type II materials through the minimization of the Ginzburg-Landau 
functional, see for instance the papers [H [H [20l [III [28l [30] and the books O [29] 
for results and additional references concerning the subject. 

In addition to the phase transitions associated with minimizers (stable states) 
of the Ginzburg-Landau functional, type II materials posses hysteresis associated 
with local minimizers (locally stable states) of the energy. For instance, a locally 
stable state that does not posses vortices will remain locally stable in increasing 
magnetic fields up to a super-heating field, and a similar phenomenon is associated 
with a sub-cooling field associated with decreasing applied fields (see [31]). It is 
therefore natural to address a similar question when dealing with the transition of 
normal states: Does local stability persists for the normal state below Hc3, and 
for the superconducting state above -ffc.,, or will there be hysteresis? That is, we 
ask whether the transition from the superconducting phase to the normal phase 
happens at a sharp critical value of the magnetic field (= Hc^ ) . 

For type II superconducting samples with smooth boundaries and surrounded 
by the vacuum, Fournais and Helffer [8] showed that the transition is indeed sharp. 
Hysteresis is excluded in [7]. The case of domains with corners and the 3- 
dimensional case [8l [26] (see also \W\ I26j) have also been studied. The reason for 
the sharp transition is essentially the monotonicity of the first eigenvalue fJ-i{B) of 
the Neumann Schrodinger operator 



for large values of the magnetic field, a property known as strong diamagnetism 
(see [8l[9]). Here F : ^ R'' is a vector field such that curlF is a constant. 

In this paper, we address the same question — the transition from superconduct- 
ing to normal phase — but for a superconductor surrounded by normal materials. It 
is well known from the superconducting proximity effect (see |13[I14)) that the pres- 
ence of a normal material exterior to a superconductor allows the superconducting 
electron Cooper pairs to fiow into the normal material in a narrow boundary layer. 
The characteristic length scale of that layer is called the 'extrapolation length'. To 
model this phenomenon, one has to consider a generalized Ginzburg-Landau theory 
where the order parameter and the magnetic potential are not only defined in the 
superconducting material but also in the normal material surrounding it. 

For a cylindrical superconducting sample of cross section C ffi^, surrounded 
by a normal material and placed in a magnetic field parallel to the cylinder axis, 
the Gibbs free energy is given by the following Ginzburg-Landau type functional 

(see my- 



(V-iBF)^ inn, fldR'^ (d = 2,3) 




(1.1) 



Here we use the notation 



Va = V - zA, 



(1.2) 



for the magnetic gradient. 
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In the above functional, C is a bounded domain such that C and fl\n 
is the cross section of the normal material^, the Ginzburg-Landau parameter k > 
is a characteristic of the superconducting material (the ratio of two characteristic 
lengths), H > is the intensity of the appHed magnetic field, a > is a character- 
istic of the normal material that depends on the temperature and its sign signifies 
that the temperature is above the critical temperature of the normal material. Fi- 
nally, TO > is the conductivity ratio of the superconducting and normal materials. 
Minimization of the functional l|l.ip will take place over finite-energy configurations 
(t/'jA) e H^{il;C) X H^{Q;M.'^). Starting from a minimizing sequence, it is quite 
standard to prove the existence of minimizers of , see [16] . 

We will always assume that 51 and are smooth, bounded and simply 
connected. 

Notice that the normal state (0, F), with F being the unique vector field satisfying 

curlF = l, divF = inn, • VF = on dn , (1.3) 

is a critical point of the functional l|l.ip . It can also be shown that this state is the 
unique normal state up to a gauge transformation (see [16]). Configurations which 
are gauge equivalent to the normal state (0, F) will be called trivial throughout the 
paper. 

Defining the set, 

Af^'^{a, m; k) ~ {H > : Gk,h has non-trivial critical points} , (1-4) 

then it is known from [16l Theorem 4.7] (see also Theorem l2.3l in the present paper) 
that the above set is bounded. 

In connection with stability and local stability of the normal state (0, F), we also 
introduce the two sets: 

A/'(a,TO; k) = {H > : Gk.h has a non-trivial minimizcr} , 

AA'°'=(a,TO;K) = {iJ > : fi^^^K, H) < 0} . 
Here /i*^^^ (k, H) is the ground state energy of the quadratic form 

Q[k, H]{cj>) = / (|V,HF0P - dx+ f (-\W^HF<I>\'' + aK^\cl)\A dx , 

Jn Jn\n V"^ / 

(1^6) 

Since Q[k, H] defines the Hessian of the functional Gk,h at the normal state (0,F), 
we see that if H & J\f^°^{a,m; k), then (0,F) is not a local minimizer Gk.h- Hence 
we obtain the following inclusion, 

Af^°''{a, to; k) C J\f{a, in; k) . 

On the other hand, the following inclusion is trivial, 

7V(a,77i; k) C 7V^''(a, m; k) . (1.8) 

One of the main results of the present paper is the following. 



(1.5) 



^In part of the existing literature on the subject (see for example [231 [4]) fl is taken to be all 
of R^. However, in the original paper of |15| . the normal material (i.e. H \ H) has been taken 
to be bounded. We take the latter point of view, in order to avoid certain technical difficulties 
which are unimportant for our present purpose. The functional analysis necessary to overcome 
those difBculties is developed in [16) . 
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Theorem 1.1. Given a > and m > I, there exists kq > such that, for all 
K > kq, the following equalities hold, 

7V'"'(a, m; k) = 7V'°''(a, m; k) = 7V(a, m; k) . 
In the literature, it is typical to introduce the following critical fields (cf. e.g. 

mm), 

ir^,^{a,m; k) = sup J^'''' {a, m;K) , i?^^ (a, to; k) = inf M+ \ 7VA*"'(a, to; k) , (1.9) 
iJc":, (ij "i; k) = sup A/'(a, to; k) , ffr^, (a. to: k) = inf \ A/"(a, to: k) . (1-10) 

!ff^3(«''^;'«) = supA/''°''(a,TO;K) , H}^l{a,m; k) =infM+ \ 7V'°''(a, to; k) , 

(1.11) 

As a corollary of Theorem we get a result concerning equality of the above 
critical fields. 

Corollary 1.2. Given a > and m > 1, there exists kq > such that, for all 
K > Kq, the following equalities hold, 

H^^^{a,m; k) = Hc3{a,m; n) = H^cli'^' "^^ ' 

and 

Hc^{a,m;K) = Hc:^{a,m; k) = Hc^^{a,m; k) . 

Remark 1.3. In [23|, the second author has established the following asymptotic 
expansion of H^n^(a, to; k) (it can also be shown to hold for H ^'^{a,m; k)): 

H}^''(a,m;K) = — (1 + o(l)) as k oo . (1.12) 

ao(a,m) 

The constant ao{a,m) in (|1.12p satisfies ^ < Oio{a,m) < 1 for to > 1 and 
ao{a,m) = 1 otherwise, and is defined via an auxiliary spectral problem (see Theo- 
rem \3.1\ below). 



Remark 1.4. Theorem \1.1\ does not cover the regime to < 1. In this specific 
regime, we have by Remark I that the nucleation field Hc^ coincides with the 
second critical field ■ This reflects one feature of superconductors surrounded by 
normal materials, that surface superconductivity can be absent (see [15]/ However, 
we need specific tools to treat this interesting case, which are beyond the scope of 
the present paper. 

We say that the transition from the superconducting to the normal phase is sharp 
if we have equality between upper and lower fields. By Corollary 11.21 it suffices to 
verify this for the 'local fields', i.e. look whether the following equality holds: 

H^n^Ja. to; k) = II^Q'^{a, to; k) ? 

One part of the present paper is devoted to this question, which Hnks to a mono- 
tonicity problem of magnetic Schrodinger operators. 

It is shown in [23] that there exists to* > 1 (see l|3.15p below for a precise 
definition of to*) such that for to > to*, the second correction term in l|1.12p is of 
order 1 and determined by the maximal value of the scalar curvature of ft. The 
precise result in this case is the following: 

Ci (a, to) 
ao{a,m) ' ao{a,mf 

(1.13) 

where the function Ci(a, •) : [to*, +oo[i~> R-|_ is defined via an auxiliary self-adjoint 
operator (see l|3.14p ) and denotes the scalar curvature of dfl. 



TT r \ 1^ Ci(a,m) , . ( 

an a, TO rw„ n mVl^ V / 



as K — > +00 
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Our final result is that all the critical fields coincide for large k under the extra 
condition that m is sufficiently large. 

Theorem 1.5. Given fl^Q, and a > 0, there exists jtiq > and if m > mg, 
there exists kq > such that if n > kq then all the six critical fields defined in 
(|1.9p . (|1.10p and l|l.lip coincide. Furthermore, their common value is the unique 
solution H = He's (a, rn; k) of the equation 

In particular, the asymptotics (|1.13p holds for all the six different definitions of the 
critical field. 

As we mentioned earlier, the essential key for establishing Theorem 11.51 is a 
monotonicity result concerning the first eigenvalue of a Schrodinger type operator. 
Actually, let Hi{B, a) be the first eigenvalue of the following operator 

- \7 BF ■ wM"^ BF + aBVa{x) inn, (1.14) 

with 

'^'"(•^) = {i mM2\0, ^"(^) = {"a^ !nM2\f^. ^^"^^^ 
Roughly speaking, we will establish that the equation in {B, a) 

admits a unique solution provided that B is large enough and a remains close to 
ao{a, m). 

In order to prove Theorem 1 1.1) we establish the following crucial estimate, valid 
for any critical point of the functional Qk,h (i-e. solution of l|2.ip below), 

a [ \ij\^dx < [ \ip\*dx. 
Jn\n Jn 

The above estimate is non-trivial and we prove it through a detailed analysis of the 
regularity of critical points of the functional l|l.ip (see Lemma [231 and Theorem l2.4l 
below). With the above estimate on the one hand, and other weak decay estimates 
established in Lemmas 12.21 and 12.61 on the other hand, we prove Theorem 12. 7^ 
which links the existence of non-trivial critical points to a spectral condition on the 
eigenvalue (|1.7p . Theorem 11.11 is then a consequence of Theorem 12.71 

The paper is organized in the following way. In Section^ we estabHsh a necessary 
and sufficient condition for the functional G^.h to admit non-trivial critical points, 
see Theorem 12.71 

In Section [3] we establish a monotonicity result for the eigenvalue fii{B,a), see 
Theorem 13.31 As a consequence, we obtain equality of local fields stated in Theo- 
rem [T31 Theorem 13.51 and Remark 13.61 

In Section HI we discuss the decay of energy minimizing order parameters, and 
we prove that they decay exponentially away from the boundary provided that the 
magnetic field is sufficiently large. 

Finally, in the appendix, we prove an improved expansion for (|1.13p in the particular 
case when is a disc domain. 

2. Proof of Theorem 11.11 

2.1. Basic estimates on solutions. We return now to the analysis of critical 
points of the functional l|l.ip . As we mentioned in the introduction, minimizers of 
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(|l.ip exist in the space H^{fl;C) x H^{fl;M.'^) and they are weak solutions of the 
associated Ginzburg-Landau equations: 

V-^curl A = {KHy^Wmlm {ip (V - iKHA)^j) , in , (2.1) 
1/ • (V - iKHA)ijj = , curl A = 1 on dU , 

where v is the unit outward normal vector of dVl. 

In order to avoid the question of regularity concerning solutions of the above 
equation, we shall invoke energy arguments and use only the weak formulation of 
(|fT|) . More predsely, {ip,A) G H\h;C) x H^{n;M.^) is a weak solution of if 
for all (f> g H'^{n;C) and a e H'^{n;M.'^), the following equalities hold, 

(w,n{x){\/ - iKHA)ii ■ (V - iKHA)(j) 

+ K\Vaix)i> + InixM^^)'^^ dx^O, (2.2) 

/ curl (A - F)(curla) - w„,{x){KH)-Hm{i' - inH A)ip) ■ adx = . (2.3) 

A standard choice of gauge permits us to minimize (|l.ip in the reduced space 
H^in^C) X H^ijh;R^) , where 

H^-^^{h;R^) = {A e H\h;R^) : div A ^ inh, i^-VA^O on 9f2} . (2.4) 

We get as immediate consequence of the Poincare Lemma (see e.g. [32l p. 16, 
Theorem 1.5]): 

Lemma 2.1. There exists a constant C > such that for all A G H^^^{fl; M^), 

1|A-F||^,(^) <C||curlA-l||^,(^). (2.5) 
Here F is the vector field introduced in il.3\} . 

Lemma 2.2. Let {ip,A) be a critical configuration of lll.l]} . i.e. a weak solution of 
(|2.ip . Then the following estimates hold : 

mL-in)<^^ (2-6) 
||(V-iKi7A)^/>||^,(f^^ < max(l, V^)K||ViL2(0) , (2.7) 

a I li'l^dx < I \^p\^dx, (2.8) 
Jn\n Jn 

||V'lli4(o) < ll^lU^(n), (2.9) 

and 

H\\cuAA - < C\mmn)MLHn) ■ (2.10) 

Here the constant C > depends only on a, m and Q,. 



Proof. The estimate l|2.6p is rather standard and is obtained in [T^. It can be 
derived using a simple energy argument as in [5], without relying on regularity 
properties of {-ip, A). 

Inserting 4> = in (|2.2p we get, 



[ {w„rix)\{S/ - iKHA)i:\^ + K^Va{x)\i;\^ + K^lnixMl^) dx = 0. (2.11) 



n 



Now l|2.7p and l|2.8p are consequences of l|2.1ip . The estimate l|2.9p is a consequence 
of HH). 
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Now, we prove (|2.10p . Up to a gauge transformation, we may assume that 
A e A). Inserting a = A — F in l|2.3p and estimating we get. 



kH |curl(A - F)fdx 
Jn 



- ll^'lli-(r^)ll(^ - z«;i/A)V.||^,(^)||(A - F)i'\\^,^^,y 

Invoking l|2.6p . (|2.7p and applying a further Cauchy-Schwarz inequahty, we get, 

kH I |curl(A-F)pda; < 
Jn 

'«max(^l,-i=^ U\\L2(n)U\\LHh)\\-^-^\\L^(hy (2-12) 



Now, by the Sobolev embedding theorem and Lemma [2?Tl we get 

||A - F||^,(^) < CsobllA - nmin) < C||curl(A - F)||^,(^). 
Thus, we can divide through by || curl A - Fj|^2(^) in (|2.12p to get (|2.1Qp . □ 

The next theorem gives the finiteness of the critical fields i/cs and is well-known 
(see [TTl [7] for superconductors in vacuum and [16] for a setting similar to ours) . 
We give an easy (spectral) proof for completeness. 

Theorem 2.3. There exists a constant C > such that if 

K > 1 and H> Ck, (2.13) 

then the only weak solution to i|2.ip is the normal state (0,F). 

Proof. We may assume — after possibly performing a gauge transformation — that 
the stationary point satisfies {ip, A) e H^{{l; C) x Hl^^{il; M^). Suppose that ip ^ 0. 
Using 1 1 I loo < 1, we have the pointwise inequality 

|(V - mi7F)V'P < 2|(V - iKi?A)^/f + 2{KRf\K - Fp. (2.14) 

Upon integration of (|2.14p and using Lemma [221 the Sobolev inclusion ^ L"' 
and Lemma [2?T| we find 

|(V - zKHF)ij\^dx < C«2||V>||2^^^^, (2.15) 

for some constant C > 0. This implies — since ip ^ hy assumption — that the 
lowest Neumann eigenvalue ^'^ {kH) of —(V — ikHF)^ in L^(f]) satisfies 

H^{kH) < Ck^. (2.16) 

However, since ft is smooth, we have (see [2) 

^l^{kH) > C'kH , y H>Ho, (2.17) 

for positive constants C and Hq independent from n. Combining (|2.17p and (|2.16p 
yields the result. □ 

In the next theorem, we show that the analogue of (|2.8p holds for critical 
points. 

Theorem 2.4. Let (?/>, A) e H^{n;C) x H^.^in;^^) be a weak solution of (2l\} . 
Then the following estimate holds, 

\ij\^dx< [ \i;fdx. 
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The estimate in Theorem l2.4l is an essential ingredient in proving Theorem l2.7l be- 
low (which links the existence of non-trivial critical points to a spectral condition) . 
In order to prove Theorem 12.41 we need to establish some regularity properties for 
solutions of l|2.ip . 

Lemma 2.5. Let ('0,A) € H'^{n;C) x Hli^{h;M.'^) be a weak solution of (2l]} . 
Then, 

(1) (^A, A) e H^{n\ dn-X)x i/2(f2^R2) . 

(2) (^,A) e C°°(f2\9f^;C) x C°°(f2 \ 9f7 ; M^). 

Let us emphasize that the complete regularity of {ip, A) up to the boundary 
of fl does not follow from standard regularity theory for elliptic PDE, and hence 
deserves to be studied independently. 

Proof of Lemma [27E[ The interior regularity of {ip, A) (statement (2) above) is ob- 
tained through a standard bootstrapping argument, see |291 Proposition 3.8]. 

We move now to the regularity. Since A G H)^^^{Q,), the equation for A in 
(|2.ip becomes, 

— AA ~ g{x) in Vt , 
coupled with the boundary conditions, 

curl A=l, u ■ A~Q , on . 

Here 

g{x) := {KH)-^Wra{,x) (ii/s (V - iKHA)ijj) e L^{9) . 

Now A S i?^(ri) follows from the W^'^-vegxAsxiiy of the curl-div system (see [l]). 

To obtain the regularity of t/j one has to be more careful, since -0 satisfies formally 
a transmission condition on dVl, 

T^i{v ■ (V - iKHA),p) = l^r;^n{^ ■ (V - ikHA)^,) . (2.18) 

Here 

T^^ : H\n) ^ L^d^l) , TS^' ■.H\R'\n)^L\dn) (2.19) 

are respectively, the 'interior' and 'exterior' trace operators, and v is the outward 
unit normal vector of dQ. 

The crucial point is now to apply a suitable gauge transformation which makes 
the transmission condition (|2.18p independent of the vector potential A. 

Let X be the solution of the following boundary value problem, 

—Ax = in , 1/ ■ Vx = ly ■ A on dil . 

Since A G H^{rt), standard regularity theory gives x € W^'P{fl) for every p > 1. 
Let us fix a choice of p > 2 such that Sobolev embedding gives W^'P C C^'" for 
some a e]0, 1[. 

Given a smooth domain K such that fl C K C K C ^, we can extend x to a 
function x G W^'P{n) such that suppx C K (see Remark on p. 257 in [6]). 
Now, defining, 

ip = i:e'^ , B = A - Vx , 
the condition l|2.18p reads formally, 

^ja(^-V^) = ^r3'=,f(z.-V^). (2.20) 
Actually, the equation for ip becomes, 

-div {wr,i{x)'\/p) = f in , 
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with 

f{x) = [2i(B • V) + f(divB) + KH\Bf] (p 

-K2(K(a;)(p + lo(a;)|^|V) ■ 

Here Wm and Va are as in l|1.15p . Moreover, the equation for f is supplemented 
with the boundary condition, 

1/ ■ \/(p ~ on d^l . 

The obtained equation for (p is of the form studied in [22^ Appendix B], where 
L^-type estimates are shown to hold for the solutions. Let us see how we will 
implement the aforementioned point. 

Lemma [221 gives \(p\ < 1. Moreover, since p G H^{il), A e H'^{Q) and B = 
A + Vx is bounded, we deduce that / e L'^{n). Applying now Theorem B.l in [22] . 
we deduce that ip e H^{n\dn). Since V = ve"'^ and x e W^^P{h) c Ci'"(f2), 

themij e H^{n\dn). □ 

Proof of Theoreml24[ Thanks to Lemma HH (?A,A) € H^{h;C) x iJ]i^(f2;M2) 
being a solution of (|2.ip . the function ?/' G C°°(r2 \ 90), and ("0, A) is a strong 
solution of ([21]) in f2 \ 917. 
Consequently, the function, 

u = = ^4; £C°°{h\dn;R). 

It is easy to verify that, 

^AIV'P Re(^(V - inHAfi) + |(V - iKiJAjV'P , ^ x<^n\dn. 
Thus, defining the function, 

f{x) = ]^W,n{x)\{V - IkHA)^{x)\\ 

we get that u = It/i^ is a strong solution in \ d^l of the equation, 

- idiv (w^Vit) + K^{Va + lnu)u + / = . (2.21) 

Here we remind the reader that the functions and Va are defined in l|1.15p . 

Let us show that u is a weak solution of (|2.2ip in fi. To that end, we need only 
verify that 

r^^^{u-^u)^^r^f,\u-^u) inL\dn), (2.22) 

where Tqq and Tq^ are the trace operators introduced in l|2.19p . and v is the 
outward unit normal vector of dO,. 

^ Notice that Vu = 2Re(V'V^) G H'^{Vt \ dn) since V e H^{n \ dQ;C) and 
n C M^, hence the trace of v ■ Vu is well defined in the usual sense. 
Now, it is easy to verify that, 

Vm 2 Re ii^Vi') = 2 Re (?A(V - iKHA)iP). 

On the other hand, since (V', A) G H'^{Q, \ d^l;C) x iJ^(f2;R^) and weak solutions 
of (|2.ip . they satisfy in particular, 

• (V - ikHAW = ^Tg^ni'^ ' (V - inHA)^) 

in L^(9f2). This shows that (|2.22p also holds in L^{dil) and consequently m = 
is a weak solution in ft of i|2.2ip . 
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Therefore, multiplying (|2.2ip by u and integrating, we get, 

^ (^^w^{x)\Vu\^ + Va{x)u^ + ln{x)u^ + f{x)u^ dx^Q, 

with u > and / > 0. This yields the estimate of Lemma [Til □ 

2.2. Weak Decay Estimate and applications. Just as in [3], we can derive the 
following weak decay estimate. 

Lemma 2.6. Assume that a > and m > 0. There exist positive constants C and 
C such that, if (ip, A) is a solution of ^2.1]) with: 

k{H~ k) > C , 

then the following estimate holds: 

a.n > " ' M 2 



max ( 1, MUn) < ^ / , dx 

< 



k{H - k) 

Proof. The last inequality is an easy consequence of (|2.6p , so we will only establish 
the first one. 

Let X G C°°(R) be a standard cut-off function such that, 

X=l in[l,oo[, x = in]-oo,l/2]. 



Define A = 1/ ii{K. — H), and 

A simple calculation yields the following localization formula: 
/ w,nix) (|(V - tKHA)xxif - IVxaI'IV'I') dx 

= Rc / (wrn{x){V -iKHA)ij ■ (V - iKHA){xliJj)) dx , 

where Wm is introduced in (jl.lSp . 

Now, we use the weak formulation of the first G-L equation in (|2.ip and get. 



/ xiil-WMrdx. (2.23) 



iwmix) {\iV - tKHA)xxM^ - \Vxxf\M^) dx + an" / xWi^? dx 

Jn\n 

in 

The next step is to give a lower bound to the first term on the left hand side of (|2.23p . 
This is done via an elementary inequality from the spectral theory of magnetic 
Schrodinger operators (see e.g. O Thm 2.9] or [3 Lemma 2.4.1]). Actually, since 
the function x\ has compact support not meeting the boundary dfl, and Wm — 1 
in ft, the following inequality holds. 



'__Wm{x)\{\' - iKHA)x\ip\ dx > kH / {curl A)\xxip\ dx . 
n Jn 

Writing curl A = 1 + (curl A — 1) then applying a Cauchy-Schwarz inequality, we 

get, 

I^Wm{x)\{V - iKHA)x\ip\'^ dx > kH Wm{x)\x\ip\'^ dx 
Jn Jn 

-KiJ||curlA - l\\LHn)\\x\^A\hin) ■ 
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Implementing the estimates (|2.10p . I|2.9p and l|2.8p we get, for some constant c > 0, 
the following lower bound, 

> Jjxx4'\^ dx - c4tP\\L2(^n)\\x\^\\h{n) ■ (2-24) 
Upon substitution in l|2.23p and a rearrangement of terms, we deduce. 



k{H - k) \x\->P\ dx + aK Ixxipl dx < cK\\ijj\\L2(^Q-)\\x\ip\\nm) 

Jn Jn\n 

+ llx'llioo(K)A-2 f WmixM^dx^K^ I Xlm^dx. 

j{dist(x,an)<A} Jn 
Implementing again a Cauchy-Schwarz inequality, we get 

k{H ~ k) I \x\il)\'^ dx + aK^ I \xx-4>\'^ dx 
Jn Jn\n 



< c1|Vlli2(o) + llx'llioo(K)A-^ / w^{xm' dx 

J {dmt(x,dn)<\} 

w / ixi-xlm^dx, 



Jn 

where the last term on the right hand side above is negative, since < xa < 1- 

Decomposing the integral / = / IxaV'I^ + / (1 ^ Xa)IV'I^i ^ind assuming 

Jn Jn Jn 

that 

k{H - k)> 2c^ , 



we get 



- max {k{H — k) , an^) l\x\ip\'^dx 
2 Jn 



n 

< ( + max ( 1, 1) ||x'llio.(R)A-2 ) / dx . 

"mj J J{Amt(x.dn-)<X) 



Recall that A = Xj ^ K{li — k) . The conditions on x and niji — k) imply that 

'{dist(a:,afJ)<A} 



max(l,^^— ) /jxAV'l'rfa;<4max(l,-)||x'||ioo«, / \^\'' dx . 



Consequently, we get. 



1 

m 



< (4max( 1,3 ) ||x'||i»(K) + l) / dx . 

I {d\st{x.,dn)<\} 



Choosing C = max |^c^, 4 max (l, llx'llioo(R') + 1^ , we get the desired bound. □ 

The next theorem gives a purely spectral criterion for the existence of non-trivial 
critical points of Ijl.ip . 

Theorem 2.7. Given a > and m > 1, there exists kq > such that, for all 
K > kq, the following two statements are equivalent: 

(1) There exists a solution (ip,-^) of 112.1]} with Wi-'W [^-Ai^f^-j ^ 0. 

(2) The parameters k and H satisfy ^'''^\k,H) < 0, where the eigenvalue 
Ii^^\k,H) is introduced in fi.?)) . 
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Proof. It is well known that the second statement implies the first one. Actually, we 
only use (iiA* , 0) , with t sufficiently small and i/'* an eigenfunction associated with 
/i*^^) (k, i/), as a test configuration for the functional l|l.ip . The resulting energy 
will be lower that that of a normal state. Hence a minimizer, which is a solution 
of (|2.ip . will be non-trivial. 

We assume now that the first statement holds and we show that the second 
statement is true provided that k is sufficiently large. 

Thanks to Theorem 12. 3^ it is sufficient to deal with applied magnetic fields sat- 
isfying H < Ck, for some constant C > 0. 

On the other hand, if c > is a sufficiently small constant, then we can show 
that 

A*^^^ (k, H) <0 for all H < ck . (2.25) 
In order to see this, notice that the variational min-max principle gives, 

(k, H) < ^^(KiJ; n)- , 

where fi^ {kH;Q) is the lowest Dirichlet eigenvalue of —(V — inHY)^ in L^{il). 
Standard estimates on Dirichlet realizations of magnetic operators (see [18]) yield 
the existence of a constant C > such that 

Ai^(B;17) < C'max(S,l) VB>0, 

implying l|2.25p . 

Therefore, we restrict ourselves to applied magnetic fields satisfying, 

ck<H<Ck. 

Using Lemma 12.61 and a Cauchy-Schwarz inequality, we find a positive constant 
C > such that 

/ \ 1/2 

|2 ^ r<\ i IL/.II2 ^ IL/.II2 



where we use Theorem 12.41 to get the last inequality. 

Since (V', A) is a (weak) solution of l|2.ip . we get by setting = V in (|2.2p 
together with the assumption on V' that 

< AWL^isi) < - L {wm{x)\{V - iKHA)^b\^ + n'^VaixM^) dx A , (2.27) 

where Wm and Va are introduced in l|1.15p . Notice that use (|2.8p impHes that 
lL4(n)?^0. 

Implementing (|2.27p in l|2.26p . we get 



ml,^~^<C"VAK-'/\ (2.28) 



We estimate 

|(V-iKffA)V'P 



> (1 - VAk-^/^)\{W - ^KHF)i;\^ - ^-^(«if)2|(A - F)V'| = 



This yields the following estimate on A: 

A< (-^,^^\K,H)+mm(^,l-)VAK-^^^^i''{KH)]\ml,^~^ (2.29) 



\ mj y/A Jn 
where fi^ [kH) is the lowest eigenvalue of the magnetic Neumann Laplacian 



[V -iKHYy in n 
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We recall that H = 0{hi). Using the asymptotic behavior of fj.^{B) as i? — *■ (X3 
together with the estimate (|2.28p . we get: 

VA^-^^'^l^{KH)U\\l,^^^^ < C^'^^'A. (2.30) 

On the other hand, using a Cauchy-Schwarz inequality and the estimate l|2.27p . we 
get: 

r ^/2A 
{nHr JJ{A- F)^p dx < {kH)'\\A Fjl^,^-^^ . 

By Sobolev embedding and Lemma [2?H 

{.HnA F\\l^~^ < C(.i/)2||curl A - < C.^H^H^.^^) • 



Using again Theorem 12.41 we deduce that, 

(KHnA - F\\l,^~^ < Ck^MIU^^^ < CA . 

Therefore, implementing all the above estimates in l|2.29p . we get 

Knowing that A > (see (|2.27p ). we deduce for k sufficiently large the desired 
inequality, ^(^^(k, i?) < 0. □ 



Proof of Theorem \l.l[ We can now finish the proof of Theorem ll.ll as follows. The- 
orem [TT] gives Af^'^{a,m; k) = Af^°^{a,m; k) for k sufficiently large. On the other 
hand, we have the trivial inclusions J\f{a, m; k) C Af^'^{a, m; n) and Af^°'^{a, m; k) C 
Af{a,m;K). □ 

3. MONOTONICITY OF THE FIRST EIGENVALUE 

3.1. The constant ao{a, m). We recall in this section the definition of the constant 
Q!o(a, m) introduced in |23], together with the main properties of a family of ordinary 
differential operators. 

Consider the space B^(R) = H''{m) n L^iM; \t\''dt), k e N. Given a, m, a > 
and ^ G M, let us define the quadratic form : 

B^{R) B q[a,m,a;^]{u), (3.1) 

where : 

q[a,m,a;^]{u)= [ {\u' {t)\' + \{t ^ OHt)]^ ^ a\u{t)\^) dt (3.2) 

- [\u'{t)\^ + \{t~ e)w(t)pl + aa\u{t)\A dt. 
m J 

We denote by i/[a, m, a;^] the self-adjoint operator associated to the closed sym- 
metric quadratic form (|3.ip . The domain of H[a,m,a;^] is defined by : 



(3.3) 



D{H[a,m,a;^])^l^ueB\R); u,^^ G B2(K±), m'(0+) = ^m'(0_) 
and for u G D{H[a,m, a; ^]), we have, 

{H[a,m,a;^]u){t)^ { (3.4) 
I [{ii{'d! + it-0'}+aa)u\{t); iit<0. 

We denote by fii{a,m,a;£,) the first eigenvalue of H[a,m,a;£^] which is given by 
the min-max principle, 

^i(a,TO,a;4)= int r^— . (3.5) 
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We summarize in the next theorem the main results obtained in [23] concerning 
the above family of operators. 

Theorem 3.1. The following assertions hold. 

(1) Given a > and m > 0, there exists a unique ao(a, m) > such that 

inf /ii (a, TO, ao (a, to) ; ^) = . 

Moreover, aQ{a,m) = 1 if m < 1, and there exists a universal constant 
00 €]0, 1[ such that Oq < ao{a,m) < 1 ifm>l, and 

lim ao(a, m) = 8o . 

m — ►oo 

(2) If < a < Q!o(a, m), then inf ^i{a,in,a;^) > 0. 

(3) Given a > 0, there exists toq > 1 such that for all to > toq the function 
M 9 ^ 1-^ ^1 (a, TO, Q;o(a, w); ^) admits a unique non- degenerate minimum. 

(4) For all a > and m > mo, there exists eo{m) > such that if 

a e [an{a, to) — eo("^), "0(0, "^) + eo("^)] , 

iften i/je function M 9 ^ 1— > /^i(a, m, a; ^) admits a unique non-degenerate 
minimum, denoted by £,{a,m, a) . 

(5) For all m > and a > 0, inf //i(a, m, a; ^) — > os a — > ao(a, ?7i). 

Let us mention that the universal constant 60 is the infimum of the spectrum 
of the Neumann Schrodinger operator with unit magnetic field in R x , and it 
holds that ^ < 60 < 1 . We point also that the existence of the constant toq > 1 in 
Theorem l3.1l is nontrivial, and is due to a fine asymptotic analysis of the eigenvalue 
([3?5|) as TO ^ 00, see O Section 3.4]. 
For further use, we introduce the constant 

/3(a, TO, a) = inf /ii(a, m, a; ^) , (3-6) 

and the set 

M(a, TO, a) = € M : /ii(a, m, a; ^) = /?(a, to, a)} . (3.7) 

We notice as a result of Theorem l3.1l that /3(a, to, a) = if and only if a = ao{a, to), 
and in this case M(a, m, a) = {£,{a, to, a)} when m > toq. 

Let us also notice that it results from a simple application of the min-max principle 
together with well known results concerning the harmonic oscillator in R+ , 

) + a > min ^60, ^ + + 1)"^ ' V a > . (3.8) 

Furthermore, let /^'™ be the normalized eigenfunction associated with the eigen- 
value (|3.5p . We introduce the constants 

Ci(a,TO,a;0 - / - O'l/^T (OP^t + - / {t ~ mfl'T it)\^dt (3.9) 

1 ^ 1 



-2il--JI4T(0)r, (3-10) 



6i(a,TO,a;0 = / X^^^^dt^a X-^^'dt. (3.11) 

If m > Too, a fills the hypotheses of assertion (4) in Theorem [HIT] and ^ = ^(a, to, a), 
then we write simply: 

Ci(a, to, a) = Ci (a, to, a; ^(a, to, a)) . (3-12) 
61 (a, TO, a) = 61 (a, to, a; ^(a, to, a)) , (3.13) 
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and if further, a = = ao(a,m), we introduce the constant (appearing in the 
asymptotic formula l|1.13p ) 

Ci{a,m) = -— -. 3.14) 

bi[a,m,ao) 

Finally, we point out that it is proved in [231 Section 3], 

lim Ci (a, TO, Q!o) = — (1 + 6a8Q)CJ' , lim 6i(a, to, ao) = 1 , 

m — >OD m — >QC 

where > is a universal constant (see [23l Remark 1.4]). Thus, for large values 
of TO, the constant in l|3.12p is negative. 

Now, for m > 1 and a > 0, we call (Vm,a) the property below, 

{M(o, TO, a) = {£,{a, to, a)} and ^(a, to, a) > , 
£_{a,m,a) is a non-degenerate minimum point, 
Ci(a,TO, a)<0 and 6i(a, to, a) > . 

We introduce the constant to* > 1, 

TO* = inf{TO > 1 : V 7Ti' > TO , {Vm',a) holds for a = a{a,m')} . (3.15) 

In view of the result of Theorem 13. H we may define to* as above. We emphasize 
also that the constant to* depends on a, but we omit that from the notation for 
the sake of simplicity. It is conjectured in [23] that to* = 1 for all a > 0. 
Now, given a > and to > to*, we introduce 

e*(TO) = sup{e e]0, do/2[ : V a G [ao{a, to) — e, ao{a, to) + e] , {Vm a) holds} , 

(3.16) 

where do := Qq — ^ > 0. 

For the special case of a disc domain, one more constant will appear to be relevant 
(this is C2(a,TO, a) introduced below). Given a > and m > m^, the lowest 
eigenvalue of the operator l|3.4p for ^ = ^{a,m,a) is /3{a,m,a) introduced in l|3.6p . 
Hence, the regularized resolvent, i.e. the operator (^ = ^(a,TO, a)), 

I , otherwise , 

is bounded in L^(]R). Letting / = /^'™ for ^ — £_{a, to, a), Wm{t) = 1 if i > and 
yJm{t) = ^ ii t < 0, then it is proved in [23l Proposition 3.6] that the functions / 
and limf are orthogonal in L^(R), hence the integral 



/2(a,TO, a) 



/ {t-0.fMt-0.fdt + — {t-OfMt-Ofdt 



is positive. Notice that we write simply Rq for the operator l|3.17p . Now, we 
introduce the constant, 

/•oo 1 /"O 

C2(a,TO,a)=/ \f{t)\^dt + — \f{t)\'^dt~4l2ia,m,a). (3.18) 

Jo i-oo 

Since to > 1 and /jj |/p dt ~ 1, it is clear that 6*2(0, to, a) < 1. Recalling that 

ao{a,m) > 60 > ^ , 
we get for a G [ao(a, "i) — e* (m), ao{a, to) + e* (to)] that 

Q;o(a,TO) - ic2(a, TO,a) > y , do := ©o - ^ • (3.19) 
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3.2. Notation, hypotheses and announcement of Main Result. We assume 
that C M^, f2 C are open, bounded and have smooth boundaries such that 

Hen. (3.20) 

Given B > and a > 0, we denote by P[B,a] the self-adjoint operator in L^(r2) 
generated by the quadratic form 

Hl,p{n) 3q}^Q[B,am - / {wrnix)\iV - iBF)(b\^ + aBVaixM^) dx , 

Jn 

(3.21) 

where Wm and Va are introduced in (|1.15p . and 

Hj^pin) = {ue L^{n) : (V-iBF)ueL2(E2)}. (3.22) 

Notice that since ft is bounded, Hgp{il) = H^{il) and hence the form domain of 
Q is independent of B. 

We introduce further the lowest eigenvalue of the operator P[B, a]: 

MB,a)= inf ^M^. (3.23) 

By taking B = kH and a = k/H, we see the connection with the critical fields 
introduced in (|1.5p : 

fi^^\K,H) = ni{B,a). (3.24) 

Furthermore, we put, 

Xi{B,a)=ni{B,a)+aB. (3.25) 

As an application of standard analytic perturbation theory, we have the following 
proposition. 

Proposition 3.2. Given a > 0, m > and k > 0, the following one sided deriva- 
tives 

Oh Ai± [n, H) = hm , 

a \ lu \ V Ai(g + £,a)-Ai(^,a) 

Ob Xi^±{B,a) ^ hm , 

e— ►o± e 

exist for all H > 0, B > and a > , and 

dBXi,+ {B,a) < dB\i.-{B,a). 

Moreover, there exist -normalized ground states ipH+ and ipn^ associated with 
fi^^^K,H) such that, 

dHfi'±\t^,H)^K\dBXi,±{B,a) [{Va{x) + 1)\^hJx)\^ dx 

\ a=K./H / 

(3.26) 

Here Va is introduced in ^.15]) . 

Proof. Let us explain briefiy why the one sided derivatives (?BAi_±(i3, a) above 
exist; the one sided derivatives in H exist for exactly the same reason. 

As we already mentioned, for real B, the operator P[B, a] is self adjoint, has 
compact resolvent and its form domain, H^{Q), is independent from B. Then, one 
can show that there exists zq sufficiently small such that the operator 

P[z, a] = V^F • w^V^F + azVa 
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is analytic of type (B) irH Dq = {z = x + iy : a; > 0, \y\ < zq}, see [25} p. 392] for 
the definition of type (B) operators. In particular, for a given B > 0, the eigenvalue 
^i[B,a) has finite multiplicity n. Now, by analytic perturbation theory (see [ISl 
Theorem 4.2, p. 395]) applied to the family of operators 0^3 z ^ P[z, a], we get 
the existence of e > and 2n analytic functions 

{B-€,B + e)Bz^(l)j{z)e H\h) \ {0}, {B - e, B + e) 3 z ^-^ Ej (z) e M 

such that 

P[z,a] i^.iz)) = E,{z)4'j{z), E,{B) = ^l,{B,a). 
By choosing e > sufficiently small, we get the existence of j± e {1, 2, • • • , n} such 
that 

Ejj^ (z) = min Ej{z) iox B < z < B + e, 

Ej_ (z) = min Ej{z) for B - e < z < B . 
ie{i.2,--- ,n} 

With this choice, it is clear that dB^i,±{B, a) = E'j^ {B) + a. 

Furthermore, we can choose the eigenfunctions 4>jj_ (z) to be L^-normalized. Now, 
the equality l|3.26p is obtained through differentiation of the relation (in z = 0± ) 

Ai (B + z,a) ^ Q[B + z, a]((/)j± {B + z)) + a{B + z) , 

and the application of the chain rule. □ 

In this section, we shall work under the following hypothesis on the constant m: 

m> m^, (3.27) 

where m* > 1 is the constant introduced in p.l5p . 
The next hypothesis is on the constant a: 

— e*(m) < a — Q!o(a, m) < e*(m) , (3.28) 

where the constants Q!o(a, rn) e]0o, 1[ and e*(m) > are introduced in Theorem l3.ll 
and l|3.16p respectively. 

Theorem 3.3. (General domains) 

Under the hypotheses p.27p and l|3.28p . if fl is not a disc, the following holds 



B 



lim sup |9i3Ai,±(B,a)-/3(a,TO,a)~a| =0. (3.29) 



■+00 



|q — ao(a,m)|<e*(m) 



Theorem 3.4. (Disc domains) 

Assume that Q, = 15(0,1) is a disc. Given a > and m > to*, there exists a 
constant Bq > I and a function [Bq,co[3 B t-^ g{B) satisfying lim g(B) ~ 0, 

such that if B > Bo and a satisfies 113. 28\} . then 

9sAi,+ (B,a) >a - ^C2{a,m,a) + g{B) . (3.30) 



In view of l|3.26p , we get as corollary of Theorems 13.31 and [ 

Theorem 3.5. Let a > and to > m*. Assume that C has a smooth and 
compact boundary, and if fl is a disc, assume in addition that, 

- ( aaoia,m) + ^C2{a,m,ao{a,m))) (l ^ — -) (3.31) 



Q;o(a, m) 

< ao{a, to) — -6*2(0, TO, ao{a, m)) . 



^Actually for zq small, P[2,«] is sectorial for all z £ Do, its form domain is independent from 
2 and the expression of the quadratic form associated to P[z,ci] is analytic in z. 
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Then there exists kq > such that, for all k > kq, there is a unique H^{k) > 
solving the equation 

= 0, 

where fi^^\K,H) is introduced in (L7^. 

Moreover, for all H > H^{k) and k > kq, H) > 0. 

Remark 3.6. Thanks to the asymptotic result of Theorem \3.1\ the condition iS.Sl]) 
is fulfilled for large values of m. Thus Theorem \3 . 1\ and Theorem \3.5\ imply Theo- 
rem \1.5[ 



Proof of Theorem \3.5[ From [23], we get a constant kq and a positive function S{k) 
such that Hm„_,oo S{k) = and 



< kS{k) . 



VKe[Ko,Oo[, fl^^\K, H) = =^ H ^ 

ao[a, m) 

It is moreover proved that one may find a solution H > such that fi^^^K, H) = 0. 

Let (k) = mill {H > Q : = 0}. It is sufficient to show that 

/i(^) (k, H) > for any magnetic field H satisfying 

(k)<H< — ^ — - + kS(k) . (3.32) 
ao[a, m) 

Thus, for such a magnetic field, we put B ~ B{k; H) ~ kH and a — a{K; H) = k/H 
so that we can pick kq > sufficiently large such that the hypothesis l|3.28p is valid, 
and hence the results of Theorems 13.31 and 13.41 hold. 

Let us show that the function H /xi(k, H) is strictly increasing in the interval 



k5{k) . 



ao{a,m) ' ao(a,m) 

We assume that k> kq. From l|3.26p . we write 



k(5(k) 



(3.33) 



dHti^±\K,H)^K(^dBXiMB,a)-aJjVa{x) + l)\cf>B^a\''dx^ , (3. 



34) 



where 4>B,a is an L^-normalized eigenfunction associated with the lowest eigenvalue 
fii{B,a). 

Assume first that fl is not a disc. Invoking Theorem 13.31 we get, provided that kq 
is large enough. 



dHfi^±H'^,H)>Ka(^l- JjVa{x) + l)\cl)B,a\^dx + o{l)^ , (3. 



35) 



where we use also that for H in the interval l|3.33p . I3{a,m,a) = o(l) as k ^ oo 
(see Theorem [XT]) . 

Since 4>B,a is normalized in L^(ri) and an eigenfunction associated with Hi{B, a), 
we write, 

a[l~ JjVa{x) + l)\(bBM^dx^ 

= a \(j)B.a\'^ dx — aa \(f>B.a\'^ dx 
Jn Jn\n 

= B-^ (^J_^w„,{x)\{V - iBF)c^B,a\^ dx ^ fii{B,a)^ . 
Using the min-max variational principle, we infer from the above. 



a 1 



UVa{x) + l)\cf>B.a\^dx) >S-i fl/.^(B)+Aii(B,a)' 
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where fi'^{B) is the lowest eigenvalue of the Neumann magnetic Laplacian 
-(V-iBF)2 in Q. 

It is a standard result that ^i'^ {B) = OqS + o{B) as B ^ oo (see [H]), and un- 
der the condition that H remains in the interval p.33p . it is proved in [23] that 
IJ,i{B,a) = o{B) as i? cx) (see Proposition 13 . 71 below) . Thus, provided that H is 
in the interval p.33p and kq is large enough, we get 

« L{Va{x) + lMB,c.\^dx] > ^+o(l), 



771 

and upon substitution in p.35p . we get when Q is not a disc 

^H^l'i\K,H)>0, Vk>ko. 

Now, we assume that is a disc. In this case. Theorem 13.41 and l|3.34p together 
yield, 

dHH±HK,H)>Kia-l-C2{a,m,a)-a{a + l)[ \(j)B,a\'^ dx + o{l)] . (3.36) 
y 2 Jn\n ' J 

Again, using \(f>B.a\'^dx = 1, we write, 

a — -C2(a, 777, a) — a(a + 1) / \(j)B,a\'^dx 
2 Jn\n 

a- 77^2(0,777, a)) I \(j)B,af dx - ( aa + -C2ia,m,a)] [ \(j)B,a\'^ dx . 



^ J .la \ ^ 



Now, we use that (f)B,a is an eigenfunction associated with fii{B,a). From the 
identity Q[B, a]{(/)B.a) = l^ii{B,a) and the min-max principle, we deduce, 

B,c.\'dx<(l-i^^^^^) [ 

n\n \ ^ / Jn 

Here fj,^{B,il) is the first eigenvalue of the Neumann magnetic Laplacian — (V — 
iBA)^ in ft. Using again the asymptotic behavior of n'^ {B) as i? 00 and the 
condition that H remains in p.33p . we deduce. 



Let us define. 



. |0i3.arrfa;< 1 ^^+0(1) / \cj)B,c.rdx. 

n\n \ ao(a,777) 



So ~ a — — C2(a, 777, Q;o(a, to)) 

-- (aQ;+ ic2(a, TO, ao(a, 777)) ) (1 - 



ao{a, m) ^ 

Under the condition l|3.3ip . 5q > Q. Under the hypothesis that H remains in the 
interval (|3.33p . we get provided that kq is large enough, 

1 f 5 

a - -C2{a,m,a) - a{a + 1) \4>B,a\'^ dx > , 
^ Jn\n ^ 

for all K > Kq. Substituting in l|3.36p . we get the desired result that the function 
H t-^ /i(^^(K, H) is strictly increasing in the interval p.33p for all k > kq. □ 

To prove Theorem 13.31 we deal separately with the case where the domain is a 
disc and the case where it is not. This will be the subject of the next sections, but 
we review first some of the known results concerning fJ.i{B, a). 
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3.3. Some facts concerning a). Let us recall the main results obtained in 

[23] for ^i{B,a), the bottom of the spectrum of the operator P[B,a] associated 
with the quadratic form p.2ip . We mention that in [23] we deal with ^^^\n,H) 
rather than ^i(B,a), but due to l|3.24p we get equivalent results for ^i{B,a). 
We start by giving the leading order term of a) as S — > oo ([23l Section 5]). 

Proposition 3.7. Given a > and m > 0, there exist a constant Bq > and a 
function [Bo,(X}[3 B t-^ 9{B) G M+ satisfying lim g{B) ~ such that, for any 

" B—^oo 

a > fulfilling f3.28\} . the following asymptotic holds: 

a) ~ (i{a, m, a)B\ < g{B)B , V B > , 

where (3{a,m,a) is introduced in f3.6\} . 

Letting (l)B,a be a ground state of P[B, a], it is proved in [23l Section 6] that for 
m > 1, 4>B,a is localized near the interior boundary 951 in the sense precised in the 
next lemma. 

Lemma 3.8. Under the hypotheses of Proposition [37% given a > 0, to > 1 and an 

integer fc G N, there exist positive constants B^ and Ck such that, 



/[dist(x,c)r!)]^|0B.a(x)pda;< , (3.37) 

Jn ^/B 

/[dist(x,9r!)]^|(V-zBF)(/.B,,(a;)pdx< — %^||0s,a|l'2(r„ , (3.38) 
Jn a/r 



for all B>Bk. 

Under the additional hypotheses l|3.27p and (|3.28p . i.e. a close to the constant 
Q;o(a,TO) and to strictly larger than the constant to* > 1 introduced in p.l5p . we 
get a two term asymptotic expansion for the eigenvalue ^i(-B, a) (see p3l Proposi- 
tions 7.1 and 7.3]). 

Theorem 3.9. Assume that a > and to satisfies f3.27]) . There exist a constant 
Bq > and a function [Bq, co[i-^ niB) G satisfying lim g(B) = such that, if 

B—^oo 

a fulfills f3. 28]) . then the following asymptotic holds: 

Mi(B,a)- (/3(a,TO,a)B + Ci(a,TO,a)(Atr)max\/B)| < g{B)VB , V B > Bq , 

where l3{a, to, a) is introduced in i3. 6\} . Ci(a, to, a) the negative constant introduced 
in f3.12\) . Kr denotes the scalar curvature of d^l and (/vr)max denotes the maximum 

of Kr. 

It should be noticed that there exist negative constants ci and C2 such that, for 
a satisfying l|3.28p . 

ci < Ci(a, TO, a) < C2 < . 

3.4. Proof of Theorem l3.3t General domains. We assume in this section that 
the domain fl is not a disc, still with smooth and compact boundary. Thus, the set 

Il = {xedn : K,{x) ^ {K,)n,^^} (3.39) 

is not identical to dfl, i.e. 11 ^ dfl. 

In this case, we can deduce from the asymptotic expansion in Theorem 13.91 that 
any ground state 4'B,a of P[B, a] is localized near the set 11. One weak version of 
expressing this rough statement is through the estimate of the following lemma, see 
[T8t [20] for sharper statements. 
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Lemma 3.10. Under the hypotheses of Theorem \ 3. 9\ for all eo > and N ^ N, 
there exist constants Bq > and Cn such that, 



J{dist(x,n)>co} 



Sketch of Proof of Lemma \3.1(A We provide some details concerning the derivation 
of the estimate of Lemma 13.101 Using the decay estimates of Lemma 13. 8[ one can 
get the lower boundd, 



Q[S,a](0s,a) > / Ub{x)\(^bAx)\ dx, (3.40) 
Jn 

where 

, , J (7-a)B, if dist(a;,9r2) > 

Ub{x) - I ^ CiK,{p{x))B^I^ - CaB^'^ , if dist(a;,9ri) < 2^-1/6 . 

Here 7 e]0, 1[ is a given constant, f3 = (3{a,m,a), Ci = Ci{a,m,a), the constant 
Co > depends on a > 0, > to* and 7, and p{x) e dil is defined by \p{x) — x\ = 
dist(x, dVL). 

On the other hand, it follows from Theorem 13.91 that 

Q[S,a]((/)B,„) < 

f/3(a,TO,a)B + Ci(a,TO,a)(Kr)maxB'/'+.g(5)Si/2A \^j^^^\^dx. (3.41) 



Combining l|3.40p and l|3.4ip . one can prove the following exponential decay using 
standard (Agmon-type) techniques, (see [24^ Theorem 6.5.4] for details) 

exp (^5oS^/''dist(a;,n)^ (I)bm < ^^||0B,a|lL2(n) 7 

where (5o > and M are constants depending on a and m. Now the estimate of the 
lemma is just a consequence of the above exponential decay of 4'B,a- D 

We start now by following ideas of Fournais-Helffer in [9] to prove monotonicity 
of the eigenvalue 

Xi{B,a)^Hi{B,a)+aB. (3.42) 
Notice that Xi{B, a) is the eigenvalue of the operator 

P[B,a] = P[B,a]+aB 
associated with the quadratic form 

<j> ^ Q[B, am = Q[B, am + .(^) , (3.43) 

where Q[B,a] is the quadratic form introduced in l|3.2ip . and both operators 
P[B, a] and P[B, a] admit the same ground states. With this point of view, it 
is more convenient to adapt the proof of |^ . 
We proceed to prove Theorem 13.31 



^ This is a consequence of combining Proposition 5.3.3 and Theorem 4.3.8 in |24) . More details 
are given also in |241 pp. 179-180] 
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3.4.1. Gauge transformation. Using an idea from we will apply a gauge trans- 
formation and work with a new magnetic potential A (instead of F) . In the new 
gauge, we will have that |A(/)b_q,| is small in the L^-sense as i? ^ oo. Notice that 
by Lemma [3.1 01 o!)R ^ is localized near the set 11, so it suffices to find a gauge where 
A = on n. 

To this end, we introduce adapted coordinates near the boundary of fl. For a 
sufficiently small to > 0, we introduce the open set 

n{to) = {x G R2 : dist(a;, dn) < to}. 

Let s 1-^ 7(s) be the arc-length parametrization of dfl and i^{s) the unit inward 
normal of dO, at 7(s). 

When to is sufficiently small, the transformation 

x] - to, to[3 (s, t) ^ 7(s) + tiy{s) e n{to) (3.44) 

becomes a difeomorphism whose Jacobian is \D^\ = 1 — tKi-{s). For x G ri(io), we 
put 

$-i(x) = {s{x),t{x)) 

and we get in particular that 

t{x) = dist{x,dil) in n, t{x) = ~dist{x,dfl) outside O . 
.Just as in [9l Lemma 2.5], we get the following lemma. 
Lemma 3.11. Let < e < \ min(<o, \dO.\), xo = $(so, 0) G and 

f7(e. So) = {x = $(s, i) G f2 : \t\ < e, |s-So|>e}- 

Then there exists a function Lp G C^(R^) such that A = F + V(/7 satisfies 

\A{x)\<CAi&i{x,dn) in r2(e, sq) , 

where C > depends only on Q,. 
Proof of Theorem\SM 

We have now all the prerequisites needed to apply the argument of Fournais-Helffer 
[9]. We include the details for the reader's convenience. 

Recall that, for e > sufficiently small and z G [B, B + e[, we associate to Ai(z, a) 
an analytic branch of eigenfunctions z i— > (ji^.a such that 

P[2;,a]0^,„ = Ai(2:,a)(/)^,„ V z G [S,B + e[. 
We may also in addition assume that |i(/)2,Q||^2^j:j^ = 1 for all z G [B,B + e[. 

Since is not a disc, the set 11 in (|3.39p is different from fi, and we can find 
eo G ^ min(io, \ d^\) and sq € such that 

[so-2eo,so + 2eo]nn==0. (3.45) 

Let A and Lp be respectively the vector field and the real-valued function defined 
in Lemma [3. Ill Let P[z, a] be the self adjoint operator associated to the quadratic 
form 

H^{n) 3u^ I (u;m(a;)|(V - izA)u\^ + az{Va{x) + l)|up) dx 
Jn 

where Wm and Va are introduced in (|1.15p . Notice that the operators P[z, a] and 
P[z, a] are unitary equivalent : P[z, a]u = e*^'^P[z, a]e~'''^'''u. Hence, 

[B,B + e[3 z^^,^^ :=e*^'^0,,„ 

is an analytic branch of eigenfunctions associated to the eigenvalue Ai(z, a) of the 
operator P[z, a]. 



\dn\ \dfl\ 
2 ' 2 



TRANSITION TO NORMAL PHASE 



23 



Calculating, 

d 



2=04 



2 w„(a:;)Rc (^-iA(t>B,a, (V - iBA)4>B.a^ + a iVa{x) + a) \(^B,aix)\'^ dx 

Pirn .(d4>B+z.a 

Q[B,a\\^ d'z 



2Re 



3+ 



where the last term on the r.h.s. above vanishes, since 4>z,a is normalized with 
respect to the norm. Thus, for an arbitrary ^ G K \ {0}, we may express 
dB^i.+ iB) in the following way 

C Jn 
and invoking the min-max principle we get further when C > 

ObMAB) > M^±C:4^M^ _c / \A4>B..\'dx. (3.46) 
By Lemma [3. m we may write, 

fjA<j>B.a\^dx<C f\dist{x,dn)\^\(j>B.a\^dx+\\A\\l^^ f \<l>B.a\^dx. 

Jn Jn ^ ' Jn\n{so,ea) 

By our choice of sq and eq, ^ \ ^(so, co) is away from boundary points of maximal 
curvature. Thus, invoking Lemmas 13.371 and 13.101 we obtain a constant i?o > 
depending only on a and m such that 



(\A(j)B.a?dx<CB-'^ iorB>Bo, 
Jn 



and we implement this last estimate in p.46p . Now, choosing C. = rjB in l|3.46p 
with 77 e]0, 1[ being arbitrary, we get from Proposition 13.71 

dBXi,+ {B) > /3(a, TO, a) + a (3.47) 

-^g((l + ^)B)--giB)^Cv, 
V V 

where g is independent from a, and g{B) as B 00. 

Applying the same argument to the left derivative dB^i.-iB), we get (the in- 
equality being reversed since z < in this case), 

aBAi-(B) < /3(a,TO,a)+a (3.48) 

+ l±^g{{l + ,j)B) + -giB) + Cv. 
77 77 

Recall that analytic perturbation theory gives dBM,+ {B) < dB^i.-iB) for all B. 
Therefore, l|3.47p and p.48p when combined together yield, 

Hmsup sup \dBXi,±{B,a) ~ P{a,m,a) — a\\ < Cr] . 

B — *C30 \ |a— ao (a.m) I <e^ (m) / 

Taking 7; ^ 0+ above, we get 



B 



lim sup \dBXi,±iB,a) — P{a,m,a) — a\ =0. 



-+00 



I Q — ao(a,m) | <e* (m) 



□ 
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3.5. Proof of Theorem 13. 4t Disc domains. In order to handle disc domains, 
we need a refined asymptotic expansion of the first eigenvalue ^i{B,a) as B — > oo. 
Let us introduce some notation. Given a > and m > to*, we introduce, 

5{n,B) =n-]^B - £,{a,m,a)\fB . (3.49) 

Theorem 3.12. Assume that Vl = D{0, 1). Given a > and m > m*, there exist 
a constant Bq > 1 and a function [Bo,oo[3 B t-^ g{B) satisfying lim g{B) ~ 0, 

and if a satisfies 113. 28]) . there exist real constants 

Sa = Soia,m,a) , Co = Co(a, to, a) , 
such that, for all B > Bq, the following expansion holds, 

fii{B, a) — ^/3(a, to, a)B — Ci(a, to, a)V^ + C2(a, m, a) (A^ + Co)^ 

< g{B) . (3.50) 

Here C2(a,TO, a) is introduced in i3.18\} and 

Ab = inf \S(n, B) — (5o(a, q;)| . 

The proof of Theorem 13.121 is very close to that of Theorem 2.5 in [8] and relies 
strongly on the fact that the eigenvalue p.23p admits a non-degenerate minimum 
in ^. For the convenience of the reader, we show in the appendix how the proof of 
[8] gives Theorem [3321 

Recall the eigenvalue Xi{B,a) introduced in l|3.42p . In view of the result of 
Theorem 13.121 above. Theorem 13.41 follows directly as [HI Theorem 2.5]. We omit 
thus the details. 



4. Exponential decay of order parameters 

The objective of this section is to prove that, for m > 1, energy minimizing 
order parameters decay exponentially away from the boundary. This estimate will 
be useful for forthcoming works dedicated to finer properties of energy minimizers. 

The main theorem is the following. 

Theorem 4.1. Assume that a > 0, m > 1 and let b > be a given constant. There 
exist positive constants M , C, e and kq such that, if (?/', A) is a solution of (|2.ip 
and the magnetic field verifies 

H 



> 1 

K 



then 



^2eV^t(x) A^|2 ^ 1 |(y _ inH A)ijA dx < C [ dx . (4.1) 

a \ J J{t{x)<^} 

Here t{x) := dist(2;, dVL). 

When TO > TO* and the magnetic field H = Hc^iK) — o(l), one should be able 
to prove a finer localization of ip, near boundary points with maximal curvature. 
Actually, an estimate similar to Lemma 13.101 for the linear problem should also 
be valid for the solution ip of the non-Hnear Ginzburg-Landau problem. However, 
to establish such an estimate will require a rather very technical work following 
previous papers |201 127j. so that we do not carry it out. Let us only notice here 
that — ^just as for the linear case — such a localisation estimate is essentially due to 
the asymptotic expansion of the first eigenvalue stated in Theorem 13.91 
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Proof of Theorem \4.1\ 

If ip = 0, then the estimate of the theorem is evidently true. Thus, thanks to 
Theorem 12.31 we may restrict ourselves to magnetic fields H satisfying 

H <Ck 

for some sufficiently large positive constant C. 

Let X S C°° (R) be a non-decreasing function with 

X=l, on [l,oo), on (-(^,1/2]. (4.2) 

Define x on O by xi^) = xW nHt{x)). Define furthermore, the weighted localisa- 
tion function / by 

fix) = x{x) exp(£v^t(a;)) (4.3) 

We calculate, using (|2.ip 

fw,n{x)(\{W - iKHA)ftlf - \Wf\''\^\^) dx + aK^ f l/V^pdx 
Jn ^ ' Jn\n 

= fm'-W') fdx 
Jn 

[ f^Hfdx (4.4) 
Jn 

Now, using Lemma Hm below, we can estimate 

Wmix)\iV - iKHA)fljj\^ dx 

= / \{W -iKHA)fij\^dx + m-'^ I \{W -iKHA)fi;\^dx 
Jn Jn\n 

>kH{1-C/\^) [ \fij\^dx. (4.5) 
Jn 

Combining l|4.4p and l|4.5p we find 

[nH-^ - C\fl^) I \f'4)\'^dx + an^ f\fipfdx < f\Vf\^\ip\^dx. (4.6) 
1 + Jn Jn Jn 

We estimate the last term 

[\Vf\^\tp\^dx<2e^KH[\f^P\'^dx + CKH[ \i;{x)\^ dx. (4.7) 

Jn Jn J {\/KHt{x}<i} 

Therefore we get, choosing e sufficiently small and for kH sufficiently large, 

i\fi^\^dx<C [ \^p{x)\^dx. (4.8) 

n J {\/KHt{x)<i} 

This implies the weighted L^-bound in (|4.ip . 



g2eVKHt|^(^)|2^^< ^ / \lb{x)\^dx. (4.9) 

n J{VKHt{x)<i} 

Inserting (|4.9p in l|4.4p (and using the same considerations) yields the weighted 
bound on (V - iKHA)t/j. □ 



Lemma 4.2. There exist constants Co, Ci such that if {ip, A) is a solution of (|2.1 
with k{H — k) > Co, then 



|(V - j«i/A)0||i.(,,) > kH{1 - , , mhin), (4.10) 

y^H{H - k) 
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for all (f) G C^{n). Also, 

for aii(j)e c^{n\n). 

Proof. We only prove l|4.10p the proof of (|4.1ip being identical. 

We estimate, using the compact support of (/) and the standard magnetic estimate 
from [2l Thm 2.9] (or [H Lemma 2.4.1]) 



||(V-mi?A)0||^2^j^2^ > Ki7 J curlA|(/.pdx 



In 

> kHUWI - {kH)\\ cm-lA - lilall^lli (4.12) 
By Lemma 12.21 and the weak decay estimate of Lemma 12.61 we have 

||curlA-l||2< f (4.13) 

H \J n{H — K) 

Furthermore, by the Sobolev inequality and scaling followed by the diamagnetic 
inequality, we find 

ml<csob{v\\yml+v-'mi) 

<Csob[v\\{^-t^^HA)^l+7r'\ml), (4.14) 

where Csob is a universal constant and 77 > is a parameter that we can choose 
freely. We make the choice 77 = 1/\/kH. Combining l|4.12p . I|4.13p and l|4.14p yields 

C . „ , C 



(4.15) 

from which (|4.10p follows. □ 
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Appendix A. Improved eigenvalue estimate for the disc 



The aim of this appendix is to prove Theorem l3.12l By assumption, Q = D{0, 1) 
and D{0, 1 + r) C ft for some r > 0. 

Let D{t) = {x € : |a;| < t} be the disc with radius t. Let Qs be the quadratic 
form 

Qb[u] = / (w„,{x)\{S/ - ^BF)^tP + aVa{x)\u\A dx , 

with domain {u e H^{D{l + r)\D{^)) \u{x) = on = ^ and |a;| = l + r}. Here 
Wm and Va are as in l|1.15p . and we emphasize that, for the sake of simplicity, we 
omit the dependence on a, m and a from the notation. 

Let Jli{B,a) be the lowest eigenvalue of the corresponding self-adjoint operator. 
Using the variational principle and the decay of the ground state (Lemma 13. 8p , we 
see that, 

Mi(5,a) =/ii(B,a)+0(B-°°) . (A.l) 
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So, it is sufficient to prove p.50p with ^i{B,a) replaced by //i(i3,a). 
By changing to boundary coordinates {s,t) (defined by l|3.44p ) the quadratic 
form Q_b[u] becomes, 



Qb[u] 

Mh 

Here 



2tt j-1/2 



,{t) {l-t)-^\{Ds-BAi)i 



\Dtu[ 



2tt f\l2 



tVa(t){\ - t)\uS^ dtds 



(1 - t)dtds (A.2) 
(A.3) 



2-K j-1/2 



{I - t)\u\^ dtds , Ai = ^-t+^ 



Performing the scaHng t 



1 , if i > 

— if t < 
m 



Vait) 



-1 , if t > 
a, if t < . 



(A.4) 



Bt and decomposing in Fourier modes, we find 
jii{B,a) = B inf es{n.B),B ■ (A. 5) 



Here the function 5{m, B) was defined in l|3.49p and eg^B is the lowest eigenvalue 
of the quadratic form qs^B on L^((— -\/Br, \/S/2); (1 — VBT)dT) (with Dirichlet 
condition, u(r) = 0, at r = —^/Br and r = ^/B/2), 



r.VB/2 
</Br 



Wr, 



)l0'(r)P 



dr 



a 



75/2 



(l-;^)K(t)|0(r)pdr. (A.6) 



Here, ^ = <^(a,TO, a) by convention (this makes sense provided that m > m*, and 
a G [ao(a,"i) — e*(m), ao(a, m) + e*(m)]). We will only consider 5 varying in a 
fixed bounded set. This is justified since it follows from [23l Proposition 4.6] that 
for all C > there exists D > such that \i\5\> D and B > D, then 

e<5,B > /3(a, ™, ") + Ci(a, m, a)B^^ + CB~^ . 

Furthermore, for S varying in a fixed bounded set, we know (from the analysis of 
the operator p.4p . especially that the minimum of l|3.5p in ^ is non-degenerate) 
that there exists a d > such that if i3 > d~^, then the spectrum of qs^B contained 
in ] — cx), (3{a, m, a) + d[ consists of exactly one simple eigenvalue. 

The self-adjoint operator [){S,B) associated to qs,B (on the space 
L^ii-VBr, VB/2); (1 - /Br)dT)) is 

M f(l - ;^)"'((^ -0+ B-HS ^ ))'] + aVa{r) . (A.7) 



We will write down an expHcit test function for [)((5, B) in (|A.12p below, giving es^B 
up to an error of order 0{B~2) (locally uniformly in S). 
We can formally develop f)((5, -B) as 



t)i6,B) = l)o + B-ni + B-'l)2+OiB—2) 



28 



S. FOURNAIS AND A. KACHMAR 



with 



dT 



+ 2(T-0(<5-^) + 2r(T-0' 



1)2 = "iWmlT) 



(A.8) 



Let uq be the known ground state eigenfunction of H[a,m, a; ^] with eigenvalue 
/3(a, TO, a). Here, by H[a, m, a; <^], we mean the operator (|3.4p with ^ = ^(a, to, a), 
considered as a self-adjoint operator on L^(M; dr). For ease of notation we will write 
l)o instead of i?[a, to, a; since they are the same formal differential operators. Let 
i?o be the regularized resolvent of f)o, which is defined by 



([)o - /3(a,TO, a)) V, J <t){T)uo{T) dr = 







<?i) II uo . 



Let Ai and A2 be given by 



:= {uo I f)iWo)L2(R.dT-) , 

'■= ^2,1 + A2,2 , 

:= {uq I f)2ito>L2(R.dT-) , 
The functions ui , U2 are given as 
ui = -i?o(f)i - Ai)iio , "2 



Ai 
A2 

A2,l 



A2,2 := ("0 I (1)1 - Ai)Mi)i2(R.rf^) 



(A.9) 
(A.IO) 



-i?o{(f)i - + (f)2 - A2)wo} 

Using the same type of argument in \21\ Proposition ILIO] or [lOl Lemma A. 5], we 
can prove that Rq preserves exponential decay at infinity, i.e. uo{t), ui{t), U2(t) 
and their derivatives decay exponentially fast as |t| 00. 
Let X G C'cr(^) be a usual cut-off function, such that 

X{t) = 1 for |t| < -g , suppx C [- 



.1 11 

4 ' 4J ' 



and let xb{t) = x(rB 4 ) . 

We define the following trial state. 



(A.ll) 



(A.12) 



iP-=Xb{uo + B ^ui+B ^^2} 
Using the exponential decay of the involved functions, we get after a calculation, 

||{[)(^,i3)-(/3(a,m,a) + AiB-^+A2i3-i)}V|L.(]_^/2,^/2[;(i-v^.)d.) (^-13) 

where the constant in O is uniform for S in bounded sets. Applying the spectral 
theorem, and noticing that (3{a, m, a) is an isolated eigenvalue for the operator [}o, 
we deduce that (uniformly for S varying in bounded sets), 

es^B = P{a, m, a) + XiB^^ + X-iB-^ + 0{B~^ . (A. 15) 

It remains to calculate Ai, A2 and, in particular, deduce their dependence on S. 
Writing, 

2r(r - 0' = {t- 0' - ie + mr - 



and using 



WmiT){T - ^)uldT = 0, 
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we get, 

Ai = Ci(a, m, a) , 

where Ci{a,m,a) is introduced in (|3.9p . In particular, Ai is independent oi 6. 

We do not need to calculate A2 explicitly. Notice that X2{S) is a quadratic 
polynomial as a function of 6. We find the coefficient to as equal to, 

/ W„(T)|wo(T)pdT-4/2, 

with 

/2 := {uo , w,„(t)(t - ORaw,n{T - (.)uq) . (A. 16) 

Therefore, there exist constants Sq^Cq G K such that 



A, = 



Recalling the definition of the constant C2{a,m,a) in l|3.18p . the above formula 
becomes, 

X2 = C2{a,m,a) {{6 - Sof + Co) . 

In Hght of (|A.ip and l|A.5p . we need only to show that 6*2(0, m, a) > to finish the 
proof Theorem 13.121 

Notice that we work under the hypothesis m > to*. This implies that ^ is the 
unique, non-degenerate minimum of the function (see p.5|) and l|3.15p ) 

z ^(z) := /xi(a, TO, a; z) . 

In particular, 

M"(e)>o. 

Now, exactly as shown in [9l Proposition A. 3], it holds that 

C2{a,m,a) = ^fi"{0: (A.17) 

yielding thus the desired property regarding the sign of 6*2(0, m, a). This finishes 
the proof of Theorem \3J2\ □ 

For the sake of the reader's convenience, we include some details concerning the 
derivation of l|A.17p . 

Sketch of the proof of (|A.17p . Let us introduce, 

E{z) = fl{z + , H{z) = --^w„,^ + Wra{t){t ^ ^ ^ zf + aVa{t) , 

together with an analytic family of eigenfunctions z v-^ (f>{z) e L^(M) such that 

ll<^(^)lli^(M) = 1, H{z)(b{z)=Eiz)<j)iz), 0(0) -lio. 

By differentiating the identity ||0(z)|p = 1 twice with respect to z, we get 

2Re((/.'(0), wo) =0, Re{q^"iO),uo) =-Wm\\ (A.18) 

Differentiating the relation H{z)(f>{z) = E{z)(/){z) we get since E{z) is minimal for 
z = 0, 

{H{0) - EiO)) 0'(O) = 2wmmt - e)"o • 
Since the functions Wm{t — and uq are orthogonal in i^(M), we get 

</.'(0) = 2 iH{0) - E{0)y^ (wr,-, it - Ouo) + cuo , (A.19) 
for some constant c e zM. 
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Differentiating twice the relation E{z) = ((f>{z) , H{z)(f>{z)) , we get, 

E"iO) = 2E{0)Re{(b'{0) , uo) - 8Re(</)'(0) , w„-, {t - ^uo) 

+ 2Re(</)'(0) , i/(O)0'(O)) + 2{uo , Wm uo) ■ (A.20) 

Substituting (|A.18P and (jA.lQP in (|A.20p . we get, 



i?"(0) = 2(^^ 



Wm(T)|uo(r)pdT-4/5 



with I2 introduced in (|A.16p . RecaUing the definition of the constant C2{a,m,a) 
we get the desired relation in l|A.17p . □ 
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